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The results based on 1992-95 data (Run 1) from the CDF and D0 experiments on the measure- 
ments of the W boson mass and width are presented, along with the combined results. We report a 
Tevatron collider average Mw = 80.456 ± 0.059 GeV. We also report the Tevatron collider average 
of the directly measured W boson width Fw ~ 2.115 ± 0.105 GoV. Wc describe a new joint analysis 
of the direct W mass and width measurements. Assuming the validity of the standard model, we 
combine the directly measured W boson width with the width extracted from the ratio of W and Z 
boson leptonic partial cross sections. This combined result for the Tevatron is Tw = 2.135 ± 0.050 
GeV. Finally, we use the measurements of the direct total W width and the leptonic branching ratio 
to extract the leptonic partial width r{W —^ ev) = 224 ± 13 MeV. 



We present new combined results on the W boson mass and width from the CDF and D0 experiments at the 
Fermilab Tevatron. We document the combination methodologies and summarize the results and the various sources 
of uncertainty, identifying those sources that produce correlated uncertainty between the two experiments' results. 
We also present the combination with the UA2 and LEP results. These measurements represent some of the main 
goals of the elcctroweak physics program at the Tevatron collider. 

The W boson mass and width arc important parameters in the elcctroweak gauge sector of the standard 
model(SM) [1]. The W boson, along with the Z boson and the photon, provides a unified description of the electroweak 
interaction as a gauge interaction with the symmetry group SU{2)l x C/(1). If this were an unbroken symmetry, the 
W and Z bosons would be massless. The mass of the W boson and its couplings, which determine its width, are 
therefore of substantial relevance to tests of the structure of the theory and the nature of electroweak symmetry 
breaking. 

In the SM the W boson mass is related to other parameters, and in the "on shell" scheme [2] it can be written as 



where a is the electromagnetic coupling constant and Gi^- is Fermi coupling constant measured in muon decay. The 
electroweak radiative correction Ar receives calculable contributions from loops containing the t — b quarks, the Higgs 
boson (which is the hypothetical agent of electroweak symmetry breaking), and any other hypothetical particles such 
as supcrsymmetric particles coupling to the W boson. Since the top quark mass has been measured [3,4], the t — h 
loop correction can be calculated. A precise measurement of the W boson mass therefore constrains the mass of the 
Higgs boson, which has not yet been experimentally observed. Should the Higgs boson be discovered in the future, 
the comparison between its directly measured mass and the indirect constraint will be a very interesting test of the 
SM. In the minimal supcrsymmetric extension of the standard model (MSSM), for example, loop corrections due to 
supersymmetric particles can contribute up to 250 MeV [5] to the predicted W boson mass. 
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In the SM W bosons decay leptonically: W ^ lu where / e {e,/i, r}, or hadronically: W — > q'q, where q,q' e 
{w, d, c, s, 6}. The leptonic partial width can be calculated [6] : 

T{W^ey) = ^^{l + 5) , (2) 

where the SM radiative correction 5 is calculated to be less than 0.5%. Including the QCD radiative corrections for 
the quark decay channels, the SM prediction for the leptonic branching ratio [7] is: 

B{W^eu)^{3 + 6[l + a,{Mw)/TT + Oial)])~^ . (3) 

Given the precision of these SM calculations, their comparison with the measured W boson width provides an impor- 
tant test of the SM. 

The precision of the W boson mass and width measurements from the LEP experiments (ALEPH, DELPHI, L3 
and OPAL) and the Tevatron collider experiments (CDF and D0) is similar, implying that our best knowledge comes 
from the combined results of all these experiments. The measurements are quite intricate with many inputs and 
incorporate constraints from data and physics models. In this situation a simple average of all measurements, with 
the assumption that they are completely independent, may be biased in the value or the uncertainty. 

In this paper we present systematic analyses of the W mass and width measurements published by the CDF and D0 
experiments at Fcrmilab. Following a brief description of the obscrvables in Section II, we discuss our methodology 
and calculations. In Section III, we consider the W boson mass as the parameter of interest and consider all other 
parameters needed for its measurement as external inputs, including the W width. We review the uncertainties on these 
external parameters as described in the respective CDF and D0 publications, identifying the relevant correlations. 
This information is used to construct the covariance matrix for combining the W mass measurements. 

In Section IV, we perform the same analysis for the direct measurement of the W boson width. This is again a 
one-parameter analysis, where in particular the W mass is treated as an external input. 

In Section V, we present a new methodology for treating the W boson mass and width simultaneously in a two- 
parameter analysis. This method departs from the previously published results in that no external information is 
used for either of these parameters. This has two significant advantages over the one-parameter analyses that are 
usually performed: (i) theoretical model dependence is reduced, making the results more meaningful and easier to 
interpret, and (ii) correlation with other methods of measuring the W mass and width is reduced, making subsequent 
comparisons and combinations more powerful. We also demonstrate that the joint two-parameter analysis results in 
no loss of precision, compared to the separate one-parameter analyses of the W mass and width. 

In Section VI, we review the analyses of the ratio {R) of W and Z boson cross sections in the leptonic channels, as 
published by the CDF and D0 collaborations. With some SM assumptions and measured inputs, this ratio can be 
converted into a measurement of the leptonic branching ratio of the W boson, and further into a measurement of the 
W boson width. We present an analysis of the correlated uncertainties and the external inputs used to extract the 
W width. Assumptions made in the extraction of the W width from R are compared and contrasted with the direct 
lineshape measurement in Section II. 

We conclude the paper with Section VII, discussing future implementations of the methodologies presented here. 
We suggest certain additional information that can be published by the individual collaborations regarding details of 
their analyses. We also mention those aspects where the collaborations may adopt analysis practices that are more 
consistent with each other. We hope that these comments will be useful for future efforts. 



II. TEVATRON OBSERVABLES 

We summarize here the observables described by CDF and D0 in their respective publications [8-14]. The 
directly measured W boson mass and width [8-11] correspond to the pole mass Mw and pole width Tw in the 
Breit-Wigner line shape with energy-dependent width, as defined by the differential cross section 

where Q is the center-of-mass energy of the annihilating partons. Cgg{Q) represents the partonic luminosity in 
hadron-hadron collisions 

A,-(Q) = — E/' -M^,Q')fj{QVs^,Q') , (5) 
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where i and j represent parton flavors, fij represent the respective parton distribution functions, x is the momentum 
fraction of the parton, and is the hadron-hadron center-of-mass energy. 

The W decay channels used for these measurements [8—11] are the eiy channel (by CDF and D0) and the fiv 
channel (by CDF). The W boson mass and width are extracted by analyzing the Jacobian edge and the high mass 
tail respectively of the transverse mass {rriT) distribution 

mT{l^,y) = ^2 pt{1^) pt{v) (1 - cos (</.(/±) - , (6) 

where pT and represent the transverse momentum and azimuthal angle respectively of the leptons. D0 has also 
measured the W boson mass by analyzing the Jacobian edge in the electron and neutrino pt distributions. The CDF 
result for the W boson mass is quoted using the mx fit, while the D0 result combines the mr fit and the lepton px 
fits taking the correlations into account. 

The W boson width is also extracted [12-14] from the measured ratio of partial cross sections 

^ _ aw B{W ^ ev) 
az ■ B{Z ee) 
_ Gw r(VF ^ ev) 

by using as inputs the calculated ratio of total cross sections, the measured Z ee branching ratio from LEP and 
the SM calculation of the partial width r{W ev). 

Eqn. 4 gives the differential cross section for the W Drell-Yan process. The extraction of Mw and Tw using Eqn. 4 
assumes the following: 

1. the W boson propagator can be described by the relativistic Breit-Wigner distribution in quantum field theory, 
and 

2. the production of W bosons in hadron-hadron collisions can be described by a factorizable process, with no 
additional interactions between the initial and final states. Since the leptonic final states are used for the 
measurements, there is no strong interaction between the initial state and the final state. Electroweak corrections 
are considered in the analyses. Higher twist effects, which in principle alter the effective partonic luminosity 
factor jCqq{Q), are in practice negligible for Q"^ = M'y^. 

3. backgrounds to the W Drell-Yan process from WW, WZ and tt production are small; in practice, these rare 
processes are further suppressed by analysis selection cuts and produce essentially no contamination. 

In addition to the above, the extraction of B{W ev) and from R assumes that the W and Z boson couplings 
to the leptons and light quarks are known, so that the inclusive cross section ratio awjoz can be calculated. The 
Z boson leptonic branching ratio is well- measured at LEP. Uncertainties associated with gw l<yz are discussed in 
Section VI. 

III. W BOSON MASS 

The Run 1 W boson mass measurements from CDF [8] and D0 [9] are 

Mw = 80.433 ± 0.079 GeV (CDF) , 

Mw = 80.483 ± 0.084 GeV (D0) . (8) 

We discuss the sources of uncertainty and classify them as being either uncorrelated between the two experimental 
results, or (partially or completely) correlated. 

A. Uncorrelated Uncertainties 

The measurement and analysis techniques used by both experiments rely extensively on internal calibration 
and collider data to measure detector response and constrain theoretical model inputs. The bulk of the uncertainty 
is therefore uncorrelated. We itemize the uncorrelated sources below. The following discussion also applies to the 
uncorrelated uncertainties in the direct measurement of the W boson width (see Section IV). 



(7) 
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• W statistics in the kinematic distributions used for the mass fits. 

• Detector energy response and resolution measured using resonances (Z, J/ip, T and tt*^). Model imcertainty 
from resonance line shapes is negligible. These data are used for the calibration of lepton energy response 
(calorimetry and tracking for electrons and tracking for muons). The Z data are also used for calibrating the 
calorimeter response to the hadronic activity recoiling against the vector boson. In the CDF analysis, the 
lepton response and resolution and the hadronic recoil are modelled by empirical functions whose parameters 
are constrained independently for the electron and muon channel. Therefore in the internal CDF combination 
of these measurements, uncertainties in the lepton and recoil models are uncorrelated between channels. D0 
performs independent empirical fits to their data which are uncorrelated with CDF fits. 

• Selection biases and backgrounds are unique to each experiment and are measured mostly from collider data, 
with some input from detector simulation for estimating selection bias. These uncertainties are uncorrelated 
between the CDF electron and muon channel measurements. CDF has no selection bias for electrons (in contrast 
with D0) because the selection cuts rely more heavily on tracking rather than calorimeter isolation, and because 
of a more inclusive W ^ ev triggering scheme. 

• The distribution of the transverse momentum (px) of the W boson is a model input, which each experiment 
constrains individually by fitting the Z boson pT distribution. Phenomenological models such as that of Ellis, 
Ross and Veseli [15] or that of Ladinsky and Yuan [16] are treated as empirical functions which, after folding in 
the detector response, adequately describe the observed pt{Z) distribution. The pT distribution is specified by 
model parameters along with Aqcd and the parton distribution functions (PDFs). The uncertainty is dominated 
by Z statistics, with small dependence on the PDFs and Aqcd- The latter introduces a small correlation between 
the two experiments which can be neglected at this level -'^ . A potentially correlated uncertainty in the theoretical 
relationship between the W boson and the Z boson pT spectra is assumed to be negligible. There is a small (3 
MeV) correlated component in the pt{W) uncertainty between the CDF electron and muon channel results. 

• The sources of background are Z ^ II where one of the leptons is lost, W ^ tu ^ e//i + i/Pv, and misidentified 

QCD jet events. The Z —t II background is estimated using individual detector simulations. The imcertainty 
on the W ^ tv ^ e//i + vvv background is negligible. The jet misidentification background is estimated by 
using loosely defined lepton data samples which enhances the background contribution (D0), or by selecting 
lepton candidates that fail quality cuts (CDF). While the techniques are similar in principle they differ in detail. 
CDF has also confirmed the jet misidentification background estimate using a photon conversion sample. The 
background uncertainties and cross-checks are statistics-limited and therefore independent. 

Table I shows the contributions to the uncertainty which are uncorrelated between the CDF and D0 measurements, 

taken from the respective publications [8,9] of the 1994-95 data as examples. All of these uncertainties should reduce 
in the future with more data, as the detector simulation and production/decay model is tuned with higher precision, 
and backgrounds are reduced with tighter cuts. 



TABLE I. Uncorrelated uncertainties (MeV) in the CDF [8] and D0 [9] W boson mass measurements from the 1994-95 

(Run 11)) (lata. W boson decay chaniicls uscci [c . ji) ar(> listed separately. 



Source 


CDF 


CDF e 


D0 e 


W statistics 


100 


65 


60 


Lepton scale 


85 


75 


56 


Lepton resolution 


20 


25 


19 


Pt{W) 


20 


15 


15 


Recoil model 


35 


37 


35 


Selection bias 


18 




12 


Backgrounds 


25 


5 


9 



'^Individual sources of uncertainty below about 3 MeV are typically not enumerated by the CDF and D0 experiments when 
the results are reported. 
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B. Correlated Uncertainties 

Sources of correlated uncertainty are associated with the modeUng of W production and decay, which we itemize 
below. The uncertainties are fully correlated between CDF and D0, with possibly different magnitudes. 

• The W boson kiiicniatic distribiitions used in the fits arc invariant under longitudinal boosts because they are 
derived from transverse quantities. The sensitivity to the PDFs arises because of acceptance cuts on the charged 
lepton rapidity. As the rapidity acceptance increases the sensitivity to PDFs reduces. The D0 W boson mass 

measurement includes electrons up to pseudorapidity jryj < 2.5, and the CDF measurement includes electrons 
and muons up to \rj\ < 1.0. The PDF uncertainty is correlated but different for the two measurements. 

• The Breit-Wigner line shape has an uncertainty due to the variation in the mass dependence of the partonic 
luminosity. This is a small contribution which D0 quotes separately, but CDF subsumes into the overall PDF 
uncertainty. 

• QED radiative corrections in leptonic W boson decays are evaluated by both experiments using the Berends 
and Kleiss [17] calculation. The uncertainty is evaluated by comparing to the PHOTOS [18] program and/or 
the calculation of Baur et al. [19]. The higher order QED effects have a different impact on the electron 
and muon channel measurements from CDF and the electron measurement from D0 due to differences in 
energy measurement techniques. We find that in the combined electron and muon channel result of CDF, the 
effective uncertainty due to QED radiative corrections is 11 MeV. This contribution is fully correlated with the 
corresponding uncertainty in the D0 result. 

• The W width input into the W boson mass measurement is provided differently by CDF and D0. CDF uses the 
SM prediction for Tw for the fitted value of Mw and the resulting uncertainty is negligible. D0 uses the indirect 
measurement of the W width which is extracted from the D0 measurement of the ratio a{W — > ei')/(j{Z ee). 
Since the line-shape fits performed by CDF and D0 for the W mass are sensitive to the assumed W width, we 
require that both experimental results use a consistent treatment of the uncertainty associated with the width 
input. For the purpose of combining the results, we take the 10 MeV uncertainty quoted by D0 to be the 
correlated error. 

Table II shows the correlated systematic uncertainties, taken from [8] and [9] respectively. 



TABLE II. Systematic uncertainties (MeV) from correlated sources in the W boson mass measurements [8,9]. 



Source 


CDF 


D0 


PDF & parton luminosity 


15 


7 ® 4 


Radiative Corrections 


11 


12 


Fw 


10 


10 



C. Combination of Results 

We use the Best Linear Unbiased Estimate [20] method, which is also used in [9], to construct the covariance ma- 
trix between the CDF and D0 measurements. For each source of correlated error, we construct a 2-component vector 
SiMw whose components are the individual uncertainties quoted in Table II, i.e. SiMw = (SiM^P^ , 5iM^) for the 
i*^ source of uncertainty. The contribution to the covariance matrix from each source is given by Vi = SiMw{SiMw)'^ , 
where T indicates the transpose. The various sources of error are assumed to be uncorrelated with each other, hence 
we add the individual covariance matrices Vi to obtain V = V^. This procedure gives us the off-diagonal term in 
the total covariance matrix V. The diagonal terms are obtained from the square of each measurement's total error. 
The square root of the off-diagonal covariance matrix element \/Vi2 gives the total correlated error between the CDF 
and D0 measurements of 19 MeV. The correlation coefficient, defined by Fi2/\/^ii^23, is 19^/(79 x 84) = 0.054. 
The combined W mass Myy for the set of two W mass measurements mi and their covariance matrix V is given by 

2 2 

Mw = {J2Hiomj)/{J2Hij) , (9) 
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where H =V and i,j run over the two W mass measurements being combined. The combined error is given by 

2 

a(A%) = (^i?.,)-i/2 , (10) 

and the for the combination is given by 

2 

= ^ K - Mw) Hij {ruj - Mw) ■ (11) 

Using this procedure, we obtain the combined result for the Tevatron colhder 

M^^™*™" = 80.456 ± 0.059 GeV , (12) 

with = 0.2 and probability of 66%. 

We note that the various W mass measurements from D0 are internally combined by D0 [9] using the same 
technique that we describe above. CDF combines its internal measurements [8] using a slightly different formulation, 
where the measurements are combined using only the uncorrelated errors, and then the correlated errors are added 
in quadrature. When the correlated errors are small with positive correlation coefficients, as we have here, the two 
formulations give very similar results. 

The result of Eqn. 12 is not very different from a simple average ignoring all correlations (80.456 ±0.057 GeV). This 
is due to the imccrtaintics being dominated by the uncorrelated components. As mentioned before, the uncorrelated 
sources of uncertainty will reduce with higher statistics. Therefore correlated theoretical errors such as QED radiative 
corrections may dominate in the future, in which case the error analysis we have presented here becomes more 
important. 

The combination of the Tevatron collider average with the UA2 measurement [21] of 

M^^2 ^ 80.36 ± 0.37 GeV (13) 

with a common uncertainty of 19 MeV yields 

M|;f = 80.454 ± 0.059 GeV . (14) 

Here we have taken the correlated component of the uncertainty between CDF and D0 as being fully correlated with 
the UA2 result, since all three hadron collider measurements are sensitive to the PDFs, QED radiative corrections 

and W width in much the same way. 

Further combination with the preliminary LEP average [22] of 

M^^^ = 80.412 ± 0.042 GeV (15) 

assuming no correlated uncertainty gives 

M^°'^'^ = 80.426 ± 0.034 GeV (16) 

as the preliminary world average (with = 0.34 and 56% probability). Figure 1 shows the W boson mass results, 
compared with the indirect value of 80.380 ± 0.023 GeV. The latter is obtained from a fit to all Z-pole data and the 
direct top mass measurements [22], as interpreted in the context of the SM. 



IV. W BOSON WIDTH 



The direct measurement of the W boson width is made by analyzing W boson candidate events with transverse 
mass above the Jacobian peak, which occurs for tot ~ 80 GcV. The fitting range extends roughly between 100 GeV 
and 200 GeV, where the resolution effects from the Jacobian peak are small. The W boson width analysis shares 
most of the issues of W production and decay modeling and the detector response with the W boson mass analysis, 
and the sources of uncertainty are therefore similar. 

As with the W boson mass analysis, the model parameters are constrained by analysis of internal data by each 
experiment separately. Therefore most of the uncertainties (shown in Table III for the 1994-95 data [10,11] as 
examples) are uncorrelated. These uncertainties are also uncorrelated between the CDF electron and muon channel 
results. 

The correlated sources of uncertainty are 
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Standard Model Indirect Prediction 



UA2 (1992) 

DO 
CDF 

Hadron Collider Avg 

preliminary LEP2 Avg 
preliminary World Avg 



= 80.36 ± 0.37 



= 80.483 ± 0.084 



= 80.433 ± 0.079 



= 80.454 + 0.059 



M^v = 80.41 2 ±0.042 



^ M^ = 80.426 ± 0.034 



80 80.25 80.5 80.75 

(GeV) 

FIG. 1. Direct measurements of the W boson mass compared with the SM prediction [22] based on a fit to all Z-po\e data 
and the direct top mass measurements. 

TABLE III. Uncorrelated uncertainties (MeV) in the CDF [10] and D0 [11] W boson width measurements from the 1994-95 
(Run lb) data. W boson decay channels used (e, ^) are listed separately. 



Source 


CDF fi 


CDF e 


D0 e 


W statistics 


195 


125 


142 


Lepton energy scale 


15 


20 


42 


Lepton E or pr non-linearity 


5 


60 




Recoil model 


90 


60 


59 


Pt{W) 


70 


55 


12 


Backgrounds 


50 


30 


42 


Detector modeling, lepton ID 


40 


30 


10 


Lepton resolution 


20 


10 


27 


Parton luminosity slope 






28 



• Parton distribution functions - the CDF and D0 analyses used different sets of PDFs to evaluate this uncertainty 
and quote different contributions. The W boson acceptance is similar in the direct measurements of the W boson 
width since both experiments require lepton pT > 20 GeV and \ri\ < 1. 

• W boson mass 

• QED radiative corrections 

The Run 1 direct W boson width measurements from CDF [10] and D0 [11] are 

Tvi/ = 2.05 ±0.13 GeV (CDF) , 

Tw = 2.2311°;}^^ GeV (D0) , (17) 

where the total uncertainty is quoted. The correlated uncertainties for the two measurements are shown in Table IV. 
The likelihood fit returns a slightly asymmetric statistical error for the D0 result. We symmetrize it by taking the 
arithmetic average and combine in quadrature with the total systematic uncertainty to obtain a total uncertainty of 
173 MeV for the D0 result. We use the procedure described in Section IIIC to construct the covariance matrix, and 
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TABLE IV. SystGHititic imcBrtciiiitiGS (M!gV) from corrGl&tGd 


sources in the direct W boson width measurements 


flO 111 


Source 


CDF 


D0 


PDF 


15 


27 


Radiative Corrections 


10 


10 


W boson mass 


10 


15 



use it to obtain the combined result 



pTovatron ^ 2.115 ± 0.105 GeV 



(18) 



with — 0.7 and probabihty of 40%. The square root of the off-diagonal covariance matrix element gives the total 
correlated error of 26 MeV and a correlation coefficient of 0.03. As in the case of the W mass combination, the 
uncorrelated errors dominate with the current statistics, and ignoring the correlation would produce a similar result 
(2.115 ± 0.104 GeV). However, in Run 2 at the Tevatron, which is expected to increase the statistics by a factor of 

20, the correlated uncertainties on the theoretical inputs may dominate. 

Combination of the Tevatron average with the preliminary LEP average [22] of 



= 2.150 ±0.091 GeV 



assuming no correlated uncertainty gives 



eWorld 
W 



= 2.135 ± 0.069 GeV 



(19) 
(20) 



as the preliminary world average (with — 0.063). Figure 2 shows the W boson width results, compared with the 
SM prediction of 2.0927 ± 0.0025 GeV [23]. 



Standard Model 
Prediction 



DO 

CDF 



Hadron Collider Avg 

preliminary LEP2 Avg 
preliminary World Avg 

I I I I I I I I I L 

1.6 1.8 2 



rvv = 2.231 +0.173 



rvv = 2.050 + 0.130 



r^v = 2.115 + 0.105 



rvv = 2.1 50 + 0.091 



rvv = 2.135 + 0.069 



2.2 2.4 

Tw (GeV) 



FIG. 2. Direct measurements of the W boson width compared with the SM prediction [23]. 



V. JOINT ANALYSIS OF W BOSON MASS AND WIDTH 



In this section we describe the analysis for the joint direct measurement of the W boson mass and the width. 
We do not allow external constraints on the mass and width parameters: instead we propagate the uncertainties on 
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the direct observables to the uncertainties on the extracted Brcit-Wigner parameters. This procedure will give us the 
values and uncertainties on and extracted from "Tevatron data only" , as well as their covariance. 

We introduce the following terminology to distinguish between the observables called Mw and Tw (which are 
returned by the fits to the data spectra) and the Breit-Wigner parameters of the same names (which we want to 
extract). We define the vector of observables o = {M^r, ^^r) and the vector of Breit-Wigner parameters t = {MlY,r\y). 
We approximate the functional dependence o ( t ) by a linear dependence, so that o and t are related by a linear 
transformation. For the error analysis, we are interested in transforming the variations in o to variations in t. This 
transformation is given by the matrix of derivatives A = do/dt, such that 



So = A6t . (21) 



The matrix of derivatives A is defined as 



A = 



(22) 



The values of the matrix elements of A have been published by the CDF and D0 collaborations, using their Monte 
Carlo simulation programs [9-11]. These simulation programs generate W bosons according to the calculated mass, 
rapidity and Pt distributions, generate the decay products according to calculated angular distributions, and subject 
them to parameterized detector response functions. The simulated decay leptons are used to predict the distributions 
of the observables in the data. 

The simulation and fitting programs demonstrate that the diagonal elements of A are unity. The off-diagonal 
element -gyr^ is given by ^ the 10 MeV variation in observed Mw due to a 60 MeV variation in [9]. The off- 
diagonal clement -^j^f- is given by the mean variation of 13 MeV ^ in observed Fw for a 39 MeV variation in My^r 
[10,11]. Thus A is reported by CDF and D0 to be 



1 0.17 
0.33 1 



(23) 



We invert Eqn. 21 to obtain A ^ So = 5t and take the expectation value of the product of each vector and its 
transpose 

A-i <dd{So}^> (A-y =<Stidt)'^> , (24) 

where T denotes the transpose and < ... > denotes the expectation value. The left-hand-side of Eqn. 24 contains the 
covariance matrix of the observables < 6o (So)'^ >, and we identify the right-hand-side with the covariance matrix of 
the extracted Breit-Wigner parameters. 

The diagonal elements of < So {Sof > are given by the variances of the individual Tevatron averages of the 
direct W boson mass and width (see Eqns. 12 and 18), excluding the error contribution to Mw due to Tw and vice- 
versa. In order to evaluate the off-diagonal matrix clement, we analyze the various contributions to the respective 
variances. The observables are obtained from fits to disjoint data samples'*, so that their statistical uncertainties are 
uncorrected. However, the observed values of Mw and Tw depend on the same detector parameters (such as energy 
scales and resolutions) and the same theoretical parameters (such as parton distribution functions and QED radiative 
corrections) . Hence the uncertainties in these "nuisance" parameters propagate into correlated uncertainties between 
the observables. 



^We use the value of the derivative quoted by DO, since CDF does not quote it. We assume that the same derivative would 
apply for both experiments since the W boson kinematics and experimental resolutions are similar. 

^Tho uncertainty in Fw due to Mw is quoted as 10 McV and 15 MoV by CDF and D0 respectively, which arc consistent with 
being equal given that both experiments round the quoted systematics to the nearest 5 MeV due to Monte Carlo statistics. 
The kinematics and acceptance for both experiments are very similar, hence we expect the true sensitivity to be the same, for 
which our best estimate is their average of 13 MeV. 

''The W mass fits are performed with the data satisfying jtit < 90 GeV or lepton pt < 50 GeV, while the fits for the W width 
are performed with data satisfying mr > 100 GeV. 
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To compute the off-diagonal term, we evaluate the uncertainty contribution to the observed M\y and Tw due to 
each of these nuisance parameters. The following procedure is followed: (i) remove the respective contribution from 
the CDF and D0 results separately, by setting each to zero, (ii) recompute the total error on the CDF+D0 average, 
and (iii) take the difference in quadrature between the original total error and the reduced total error. Table V shows 
the uncertainty contributions from each source to the CDF-I-D0 averages. 

In the above procedure, we have followed the same assumption that is made by CDF and D0 in their publications 
- that the sources of uncertainty listed in Table V are mutually independent. This is a valid assumption given the 
statistics of these data. The lepton energy scale and resolution are derived from the observed peak position and 
width of the Z boson mass distribution, which are essentially decoupled. The pt{W) uncertainty is dominated by the 
statistical error of the pt{Z II) measurement, although this may change in the future. The recoil model is tuned 
using transverse momentum balance in Pt{Z — > II) events, where the lepton resolution is a small effect compared to 
the recoil resolution. 



TABLE V. Correlated uncertainties (MeV) between the CDF+D0 averages of Mw and direct Tw, due to nuisance param- 
eters. 



Source 


Mw 


Tw 


Lepton scale 


37 


17 


Lepton resolution 


12 


11 


Pt{W) 


9 


24 


Recoil model 


20 


35 


Detector modeling, selection bias 


6 


13 


QED radiative correction 


11 


10 



We use the information from Table V to evaluate the covariance term 

< SM^ ST°yy> = Y^ SiM^ SiT°yy , (25) 

i 

where the sum is performed over the various sources in Table V, and 5iMy^ and SiV^ arc the respective error contri- 
butions to and F^y from source i. In this sum, the relative sign of each pair of factors 6iM^ and SiT^r determines 
the sign of the covariance contribution. The W mass and width analyses were performed by each experiment in a 
closely related manner, using the same simulation programs for both analyses. The uncertainty contributions due to 
the nuisance parameters are completely correlated between the observed Mw and Tw- Therefore 5iM°Y and SiV^ 
have the same sign in all cases. To illustrate, in the cases of the lepton energy scale, lepton energy resolution, pt(W) 
and recoil modeling, an increase in the respective parameter increases the observed values of both Mw and Tw- 
Similarly, in the cases of detector modeling, selection bias and QED radiative correction, the bias in the shape of the 
rriT or lepton pT spectrum affects both observables in the same direction. 



TABLE VI. Uucorrclatcd systematic niieertaiiities (MeV) 


between the CDF+D0 ;iv 


erag'es of Mw and direct Vw- 




Source 


Mw 




Vw 


Backgrounds 


6 




21 


PDF, parton luminosity 


12 




22 



Table VI shows the systematic error contributions due to PDFs and backgrounds. We do not expect a strong 
correlation between the error contributions to the observed mass and width from these sources, because the observables 
are derived from different ranges in mr- Thus, in the case of the PDFs, a different x range is relevant in each case. 
Furthermore, in the case of the W mass, the uncertainty in the PDFs propagates mainly through acceptance effects, 
while in the case of the W width, the main effect is through the relative normalization of the high and low ttit regions. 
In the case of backgrounds, QCD jet misidentification produces the dominant background whose shape is determined 
independently in the different regions. The sensitivity to the background shape and normalization is different in 
the fits for the mass and the width, since the shapes of the signal distributions are very different in the respective 
fitting windows. On the basis of these arguments, we take the contributions in Table VI to be micorrelated. They 
are not used directly in this joint error analysis; we present them here for completeness and future reference. 

Evaluating Eqn. 25, we find < 5M^ SV^ > = 43^ MeV^, and the covariance matrix for M^ and is 



15 



< 5o{5of > = 

(592 - 102) MeV^ 432 MeV^ 

MeV^ (1052 - 132) MeV2 



(26) 



where the removal of the 10 (13) MeV systematic on the individual measurement of (J^w) due to (M-i^) 
variation is shown explicitly. Substituting this result and A into Eqn. 24 gives the covariance matrix for the extracted 
Breit-Wigner parameters M^^r and F^^ : 

<dt{dt) >- (^_(33MeV)2 (106MeV)2 J ' ^^^> 

The negative sign of the covariance between and T\y can be understood as follows: a higher value of the 
Breit-Wigner pole mass increases the predicted number of events at high rriT, causing the inferred to reduce 
(given the number of observed events at high itit)- Similarly, a higher value of the Breit-Wigner width increases 
the expected number of events on the high side of the Jacobian edge, causing the inferred to reduce (given the 
observed position of the Jacobian edge). 

We now describe the calculation of the Mw and Tw central values in the joint analysis. We shift the observed 
value of each variable by the corresponding difference of the other variable from its assumed value, scaled by the 
appropriate partial derivative: 

M* = M° + b {Ta - r*) , 

r* = r° + a {Ma - M*) . (28) 

Ma and Ta denote the assumed values of the mass and width used in the width and mass analyses respectively. r° and 
M° are the values extracted in these individual analyses, and a and b denote the partial derivatives a = dV/dM = 0.33 
and b = dM/dV = 0.17. Solving these simultaneous linear equations for M* and F*, we obtain the central values in 
the joint analysis. 

CDF and D0 used different assumed values of the W mass and width in their respective analyses. For our 
simultaneous Mw — Fjy analysis, we need the individual Tevatron averages of Mw and Tw for which the inputs are 
quoted using the same reference values of F^^ and Mw respectively. Thus we cannot use the results of Eqn. 18 and 
Eqn. 12 for F° and M° directly. To arrive at the appropriate averages, we use the partial derivatives mentioned above 
to "shift" the CDF and D0 measurements to common reference values of Ma = 80.413 GeV and F^ = 2.080 GeV. 

This reference point is calculated as follows. In the width analysis, D0 assumed a W mass value of 80.436 GeV while 
CDF assumed a value of 80.400 GeV. We use the weights derived for combining the CDF and D0 width measurements 
(Sec. IV) to obtain the average Ma = 80.413 GeV. Similarly, in the mass analysis, D0 assumed a W width value of 
2.062 GeV and CDF assumed a value of 2.096 GeV. Using the weights derived for combining the mass measurements 
(Sec. Ill), we obtain the average F^ = 2.080 GeV. 

For these coordinates of the reference point, the CDF and D0 W mass measurements shift by about 3 MeV each, 
and the W width measurements shift by about 4 MeV and 8 MeV respectively for CDF and D0. The combination 
of these shifted values is then repeated according to the procedure described in Sec. Ill and IV. We obtain the new 
individual Tevatron averages of F" = 2.115 GeV and M° = 80.456 GeV. These values are identical to those quoted 
in Eqn. 18 and Eqn. 12, proving that our calculated reference point is consistent with the original choices made by 
CDF and D0. 

We can now solve the simultaneous linear equations given in Eqns. 28, to obtain 

^Tevatron ^ ± 0.059 GcV , 

pTevatron ^ 2.102 ±0.106 GeV , (29) 

as the Tevatron results of the joint analysis. The correlation coefficient is —0.174. 

Finally, it is of interest for future, higher precision measurements of Mw and Fi^,' to pursue this joint analysis 
technique. We expect most error contributions to scale with the statistics of the data. Assumptions that are made in 
providing external input for Tw in the Mw analysis are not necessary in this joint analysis technique. We also note 
that there is almost no loss of precision compared to the individual measurements. While this may seem surprising, 
the reason is the positive covariance induced between and F^ by the uncertainties in the nuisance parameters. 
This means that an error in any of the nuisance parameters moves Mw and Tw in the same direction. But since 
an increase in one causes the other to reduce as mentioned above, this overall negative feedback suppresses the 
systematic uncertainties from the nuisance parameters on both M^^ and T\^. This reduction in other systematic 
errors compensates for the information lost in excluding external mass and width input. 
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VI. INDIRECT W BOSON WIDTH AND LEPTONIC WIDTH 

The CDF and D0 measurements of R (see Eqn. 7) have been presented [12-14] elsewhere. We describe here 

the combination of the R measurements and the extraction of Tw from R assuming the vahdity of the SM. We also 
combine this extracted value of Tw with the directly measured Tw from the tot spectrum shape. 

A. Combination of R Measurements 

The published CDF [12] and D0 [13,14] measurements of R in the electron channel are 

R = 10.90 ± 0.32 (stat) ± 0.30 (syst) (CDF) , 

R = 10.82 ± 0.41 (stat) ± 0.36 (syst) (D0 Run la) , 

R = 10.43 ± 0.15 (stat) ± 0.23 (syst) (D0 Run lb) , (30) 

where Run la refers to the 1992-93 data and Run lb refers to the 1994-95 data. The uncertainties are summarized 
in Table VII. 

TABLE VII. Fractional uncertainties (in %) in the CDF [12] and D0 [13,14] measurements of R. "la" and "lb" refer to 
the 1992-93 and 1994-95 data respectively. The cohimn labelled "common" indicates the correlated error, taken as common 
between the D0 Run la and Run lb measurements. The la and lb columns indicate (in some cases additional) uncorrelated 
errors for the D0 measurements. The last column indicates the error components that are correlated between the CDF and 
the combined D0 measurements. 



Source 


CDF 


lb 


D0 
common 


la 


CDF & D0 
Correlated 


PDF 


1.1 




0.3 




0.3 


Mw 


0.1 




0.1 


0.7 


0.1 


Boson pt 


0.2 




0.1 


0.4 




Energy Scale 


0.4 


0.7 




0.3 




Recoil Response 


0.6 


0.2 




0.6 




Clustering Algorithm 






0.2 






Generator 






0.3 






Electroweak Corrections 


1.0 




1.0 




1.0 


Backgrounds 


1.5 


1.7 


0.1 


2.3 




Efficiencies 


1.5 


0.6 




1.9 




NLO QCD 


0.6 










Drell-Yan 


0.2 










Total Systematic 


2.8 


2.2 


1.1 


3.3 


1.0 


Statistical 


2.9 


1.4 




3.8 




Total 


4.0 


2.7 


1.1 


5.1 


1.0 



We combine the D0 results from Run la and Run lb taking the systematics due to choice of PDF (0.3%), the 
uncertainty in Mw (0.1%), the uncertainty in the boson spectrum (0.1%), clustering algorithm dependence (0.2%), 
physics generator issues (0.3%). electroweak radiative corrections (1.0%) [19] and Drell-Yan background (0.1%), as 
correlated error components, to obtain a total correlated uncertainty of 1.1%. The result for the combined D0 
measurement is 

i?D0 = 10.50 ± 0.23 (uncorrelated) ± 0.12 (correlated) 

= 10.50 ±0.26 . (31) 

This D0 result is then combined with the CDF measurement. In this combination, the systematics due to the 

choice of PDF (0.3%), the uncertainty in Mw (0.1%) and higher-order electroweak corrections (1.0%) are treated as 
correlated uncertainties, to obtain a total correlated uncertainty of 1.0%. The average R value is 



^Tevatron ^ 10.59 ± 0.20 (uncorrelatcd) ± 0.11 (correlated) 
= 10.59 ±0.23 . 



(32) 
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B. Extraction of W Boson Width 

In the extraction of Tw from R, the Z ^ ee branching ratio is taken from the PDG [7] to be (3.363 ± 0.004)%. 
The inclusive cross section ratio (Tw/uz is calculated at NNLO using Van Nccrven et al. [24], with the following 
inputs: M^?™*™" = 80.456 ± 0.059 GeV, Mz = 91.187 GeV, Tw = 2.06 ± 0.05 GeV, Tz = 2.490 GeV, and 
sin^^i^ = 0.23124 [25]. The renormalization and factorization scales are set to the boson mass. The calculated value 
of the ratio of inclusive cross sections is found to be 

— = 3.360 ± 0.051 . (33) 

The dominant uncertainties (quoted in parentheses) in the calculation of the cross section ratio are due to PDFs 
(0.45%), Mw (0.09%), factorization scale (0.12%), renormalization scale (0.06%) and sa^dw (1.43%). The uncer- 
tainties due to Mz, and Vz arc negligible. The uncertainty due to any input is estimated by varying the input 
by ±1(7 and taking half of the difference between the results. The uncertainties due to the renormalization and 
factorization scales are estimated by varying the scales high and low by a factor of two. The uncertainty due to 
electroweak corrections is estimated by taking different conventions for sin^^v^. For our central value we use sin^^g// 
from LEP, which gives an effective Born approximation and minimizes higher order corrections. The on-shell value 
of syc?Qw = 1 ^ [Mw jMzY 1 however, is equivalent at tree level, but gives a Z boson production cross section which 
is about 1.4% higher. We include this variation as a systematic uncertainty on the calculated uwjf^z- 

The value of the W ^ ev branching ratio extracted from the combined CDF and D0 measurement of R using 
Eqn. 7 is 

B{W ev) = (10.61 ± 0.28) % (Tcvatron) . (34) 
For comparison, the SM value of the branching ratio [7] is 

ViyV ev) 



B{W ev) 



T{W) 

= [3 + 6{1 + ^ + 1.409 (^)2 - 12.77 {—f}]-^ 

TT IT TT 

= 0.10820 ±0.00007 , (35) 

where we have used as{Mw) = 0.1224 ± 0.0028 ^. 

The SM calculation of the W boson leptonic partial width is given by [6] 

r(W-^ez/) = ^^(l + 

= 227.1 ±0.6 MeV , (36) 

where = (1.16637 ± 0.00001) x 10"^ GeV'^ is the muon decay constant [7], 5^^ = -0.0035 ± 0.0017 is the 
"oblique" correction to the tree level partial width [6], and M^^™*™" = 80.456 ± 0.059 GeV. Using this value of 
T{W ev), the extracted value of Tw is 

pTcvatron ^ 2.141 ± 0.057 GcV . (37) 

For comparison, the SM prediction is = 2.099 ± 0.006 GeV, using Eqns. 35 and 36. 

A comparison of the width extracted from R with the directly measured width (Eqns. 18 and 29) provides an 
interesting test of the SM, since the two methods are quite different. In the former measurement assumptions are 
made about boson couplings, whereas the latter makes use of kinematics. This test is one of the main goals of the 
Tevatron electroweak physics program in the future, as the precision of both measurements improves with more data. 

If one is willing to make all the SM assumptions mentioned in Sec. II, it is possible to combine the indirect and 
direct measurements of F^^. The result should be used with care; for instance, it may not be used in global fits where 
boson couplings are free parameters, or in analyses of data where new physics can affect the W boson branching 
ratios. 

Given these caveats, we discuss other aspects of combining the indirect and direct measurements of Tw We consider 
the correlation induced by theoretical inputs used in the respective analyses: 



^This value of as{Mw) is obtained by evolving as{Mz) = 0.1200 ± 0.0028 (Eqn. 10.50 of [7]) from Mz = 91.19 GeV to 
Mw = 80.45 GeV using Eqn. 9.4 of [7]. 
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• PDF uncertainties: we conclude that there is no significant correlation because different aspects of PDFs are 
relevant for each analysis. For the direct measurement of the width the PDFs influence the mass dependence of 
the Breit-Wigner line shape at high mass. For the extraction of the width from R, the PDF's influence the boson 
acceptance via their rapidity distributions. The total cross section ratio aw/crz is affected by the u{x)/d{x) 
ratio of PDFs. 

• Electroweak corrections, factorization and renormalization scales and siaP'Ow play no significant role in the direct 
Tw measurement. 

• We consider the correlation induced by variation in M^^. The uncertainties in R and dwl^z due to M^^ 

variation are of the same magnitude and sign^ and therefore cancel in B{W ev). The uncertainty in the 
SM calculation of T{W ev) due to uncertainty in Mw is 0.3%, which is transferred to the extracted Tw 
as a 7 McV uncertainty. This is anti-correlated^ with the corresponding M^^ uncertainty on the direct F^^ 
measurement (13 MeV). 

Taking the anti-correlation induced by Mw variation into account^, we find the Tevatron combined (direct and 
indirect) result 

pTevatron ^ 2.135 ± 0.050 GeV . (38) 

The of this combination is 0.05 with a probability of 83 %, indicating consistency between the direct and indirect 
measurements. Further combining with the preliminary LEP direct measurement (Eqn. 19) gives 

pworid ^2.139 ±0.044 GeV . (39) 

Our world average differs from the PDG [7] value of F™''''^ = 2.118 ± 0.042 GeV because we have considered the 
correlations between the CDF and D0 measurements, which were ignored in [7]. 

C. Extraction of W Leptonic Width 

We may use the extracted value of the W ^ ev branching ratio (Eqn. 34) and the directly measured total W 
width (Eqn. 18) to obtain a measurement of the W leptonic partial width 

F(W ^ eu) = Twy~ B{W ev) 

= 224 ± 13 MeV (Tevatron) . (40) 

The fractional uncertainty in the direct Tw (5.0%) dominates over the fractional uncertainty in B{W eu) (2.4%). 
This measurement of TiyV — > ev) is in good agreement with the SM calculation given in Sec. VI. 

VII. SUGGESTIONS FOR FUTURE PUBLICATIONS 

There are a few instances where CDF and D0 have treated uncertainties differently in their respective analyses. 
For future efforts it would be helpful if a consistent treatment were adopted by both collaborations with mutual 
agreement. We itemize these cases below. 

• The uncertainty due to PDFs has been treated differently in two respects, the acceptance-related uncertainty 
and the parton luminosity uncertainty. 

The acceptance-related effects were studied by varying PDFs, but these variations differed between the CDF 
and D0 analyses. CDF used their W ^ ev charge asymmetry data to constrain the PDF variation, whereas the 
variation considered by D0 did not have this constraint because D0 did not have electron charge discrimination 
capability in Run 1. With the Run 2 detector D0 can also make this measurement. D0 and CDF have also 



With increasing Mw, R reduces due to increased acceptance, and aw also reduces. 

'^With increasing Mw, the calculated r(W ev) increases, whereas the directly measured Tw decreases. 
^Ignoring the anti-correlation changes the result and uncertainty by less than 1 MeV. 
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demonstrated that the boson (decay lepton) rapidity distributions measured in Z (TV) boson events can provide 
additional PDF constraints, especially when forward lepton coverage is included. The optimal use of all this 
information would be to impose a combined constraint on PDFs using CDF and D0 data, and then propagate 
the same PDF uncertainty into their respective analyses. 

This approach may also be applied to the parton luminosity uncertainty. Furthermore we suggest that the 
possibility of correlation between the parton luminosity uncertainty and the acceptance-related PDF uncertainty 
be studied. We suggest that these components be quoted separately along with their correlation. 

• The treatment of the W width input in the W mass analysis and vice versa should be standardized. We suggest 
that both experiments adopt a common reference point based on available information. This also implies using a 
fixed width (mass) for the mass (width) fitting instead of building in a SM relationship between these parameters. 
This approach will facilitate the combination of the 1-parameter measurements, the 2-parameter joint analysis, 
and the comparison to theory. 



VIII. CONCLUSION 



We have presented the Run 1 results on the W boson mass and width from the CDF and D0 experiments, and 
examined their sources of uncertainty to identify the correlated components. We have used the covariance matrix 
technique to combine the respective measurements from the two experiments. The probability for each combination 
is good indicating that the measurements are consistent. We have also reported the values and covariance matrix of 
the W mass and direct W width measurements from their joint analysis. Finally, we have combined the measurements 
of the ratio of W and Z boson cross sections, and extracted the combined value of the W leptonic branching ratio 
and the total W width. The measurements of the W width using the direct and indirect techniques are consistent, 
providing a test of the standard model. We have also extracted the W leptonic partial width from the measured total 
W width and the leptonic branching ratio. We have documented the methodologies that can provide the basis for 
future work with data of higher precision. 
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